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Abstract 

The order - a radiative corrections to the differential decay rate of polar- 
ized orthopositronium are obtained. Their influences on the three photons 
coincidence rate as a function of positronium polarization is considered. 
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Positron experiments for QED testing are of interest at present due to discrepancy be- 
tween calculations and the most precise measurements of orthopositronium decay rate [|l],^ . 
Although the recent experiment 0] shows good agreement with the order-a theoretical 
results this problem is certainly exist and new measurements are necessary. One of the pos- 
sible kinds of experiments is the measurements of angular distribution of the decay photons 
formed by polarized positronium proposed in ref. [Q. In this experiment a three-photon 
coincidence rate will be measured as a function of decay plane orientation. It allows to de- 
termine the products of annihilation amplitudes corresponding to the different positronium 
spin states. Usual measurements of positronium decay rate only give information about the 
sum of decay amplitudes squared. 

In this Letter we have calculated radiative corrections to the differential decay rate of 
polarized orthopositronium. Until the present time the similar calculations were carried 
out for nonpolarized positronium only ||5|-[TT|. The expression for differential decay rate of 
orthopositronium is 

dT 



1 d^ki d^k2 d^h 6{P - h - k2 - h 



Arrie 2cui 2uj2 2uj^ {2tiY 
X E \Y.nr^)P^.nT\n). (1) 

Here me is the electron mass, P = (0, 2me) is the positronium rest-frame energy-momentum 
vector, ki is energy momentum vector of final state photon with energy Ui and polarization 
vector ei , pm,n is the positronium density matrix and T(m) is the decay amplitude in m 
spin state. 

When S function was integrated out expression (1) took the form 

dV = Pm,nFn,mduJiduJ2 ^'^^^'^'^ ^ (2) 
m,n 

where solid angle doi contains the first photon momentum direction, (pu is the angle of decay 
plane rotation around ki, and -Fn,m is proportional to the positronium decay matrix 

^--=2^ E jri^mm). (3) 

Polarized positronium decay is described by five parameters. Three angles define orientation 
of decay plane and energies of two photons provide information about relative orientation 
of photon momenta in this plane. Elements of decay matrix obviously depend on the choice 
of reference frame. In our calculation we fix the reference frame as follows. The z axis is 
directed perpendicularly to the decay plane. The axis x goes along the momentum of the 
first photon ki and projection of /c2 to the axis y is positive. Projection of positronium 
spin to the axis z is used to describe the positronium spin state. When choosing the label 
of photon momenta we expect the conditions ui > uj2 > 003 to be met. The choice of 
reference frame described above does not affect the generality of our results because one can 
always rewrite the positronium density matrix in our reference frame using angular moment 
transformation properties. 
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Elements of the decay matrix are not independent. Eq.(3) shows that Fij = F*. Due to 



invariance under CP-transformations [|T2[ Fi q = -F_i q = 0. The CPT invariance requires |T3 
that Fi^i = F_i_i. Thus, only two real Fi i and Fq^q and one complex number = 
|Fi _i|e"'' are necessary to describe the polarized positronium decay. 

The eigenvalues and eigenvectors of the decay matrix can be used instead of matrix 
components. This approach seems more useful because the eigenvalues and eigenvalues do 
not depend on the choice of reference frame. The eigenvalues can be expressed in terms of 
decay matrix elements Ai = Fq q, A2,3 = -Fi i ± One of the eigenvectors el = is 

always perpendicular to the decay plane and two others 

62 = exSin{ip/2) + eyCos{ip /2), 

63 = (ixC0s{ip/2) — eySin{ip/2) (4) 

lay in this plane. If one knows the eigenvalues and eigenvectors he can always represent the 
decay matrix as follows 

Fij = Aieuelj + A2e2ie*2j + AsCsiCg^-. (5) 

We have calculated the decay amplitudes as a function of photon energies and polariza- 
tions of photons and positronium. Then elements of the decay matrix were found as a sum of 
corresponding products of amplitudes over photons polarization. Similar method was used 
before by Burichenko [|T^] and Adkins |TT[] in their evaluation of order-a^ corrections to the 
positronium decay rate coming from the squares of one loop amplitudes. 

Both Feynman and three-dimensional transverse gauges of decay photons were used in 
our calculations. Virtual photon propagator was taken in Feynman gauge. The graphs con- 
tributing to the lowest order decay amplitudes (a) and order-a radiative corrections (b-f) 
are shown on Fig. |I|. Contributions of five diagrams with self-energy and vertex operators 
inside were obtained analytically. Self-energy and outer vertex graphs were evaluated us- 
ing well known expressions for one half on mass shell vertex operator and self-energy 
operator The analytical expression for off mass shell vertex was obtained to evaluate 



the inner vertex graph. Both Feynman parameters technique and direct loop integration 
in double vertex graph were used. In the latter method fco integration was done by poles 
and remaining three dimensional integrals were calculated numerically as described in 0. 
For the annihilation diagram contribution we followed the treatment proposed in |T6|,p!7 



One integration over Feynman parameters was performed analytically. The remaining two 
dimensional integrals were computed numerically. Direct loop integration was carried out to 
evaluate the binding graph. Averaging over directions specified by k and —k was performed 
for better integral convergence and terms proportional to ~ a/f were subtracted as described 
in [^,0 . Monte Carlo integration was used for numerical calculations. Gamma matrix prod- 
ucts were evaluated numerically. The program of symbolic manipulations REDUCE was 
used for computation. 

The eigenvalues of F^i in the lowest order have the form ^ 

A° = N{3R - 2C)/2, 
A° = Ar(-i? + 2C)/2, 

A? = A° + A° = NR. (6) 
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Here is the normalization constant 

N = 3To/[8ml{n^-9)], (7) 

ro = 7^(7r'-9)m«' (8) 

is the lowest order decay rate, 

R= {l-fii- n2)^ + {l-ni- n-if + (1 - ^2 ■ ^3)^ 

C = (1 - ■ n2){l - ni ■ n-i) + (1 - ni ■ n2)(l - n2 ■ n-^) 

+ (1 - ni ■ n3)(l - ^2 ■ ^3), (9) 

= ki/uJi, and their scalar products fii ■ n2 = (0^3 — uj\ — ijJ^I2uj\uj2-, etc. Angle in this 
case can be expressed as follows 

tanv3° = -2(n3^n3j,(l - nx ■ ^2)^ + n2xn2y{X - ^1 ■ ^3)^ 
- ^2 ■ ^3)^)/ 

-^(n\^-n\)(\-n2-n^f) (10) 

We assume here that photons momenta lay in the plane xy. 

We have calculated the changes of eigenvalues and the angle of eigenvectors rotation 
around the normal to the decay plane due to radiative corrections. Our results are presented 
in Table |. Their errors correspond to ±1 in the last digit. The radiative corrections cause 
significant (up to 2-5 %) decrease in the decay matrix eigenvalues in comparison with the 
lowest-order results and turn the decay matrix eigenvectors around the decay plane normal. 
Changes in eigenvalues due to radiative corrections are almost uniform in the "center" of 
phase space. The anisotropy degree rii^ = Aj/Afc changes significantly (1-2 %) near the 
"edge" of phase space only. These changes are essentially smaller than the corrections to 
the differential decay rate. Radiative corrections to the differential decay rate averaged over 
positronium polarization are in a good agreement with previous results [||. 

This work was supported by the Fundamental Research Foundation of the Republic of 
Belarus. 
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FIGURES 




FIG. 1. Feynman graphs contributing to the orthopositronium decay in the lowest order (a) 
and in the order-a radiative corrections (b-j). They are self-energy (b), inner vertex (c), outer 
vertex (d), double vertex (e), annihilation (f) and binding (g). 
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TABLES 



TABLE I. Order-a radiative corrections to the decay matrix eigenvalues and change of eigen- 
vectors rotation angle around decay plane normal due to radiative corrections. 
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0.9063 


2.37 


2 


25 


2 


39 


0.24 


0.0385 


0.9854 


4.2 


3.4 


5.6 


0.65 


0.3469 


0.8687 


2.35 


2 


13 


2 


47 


0.31 


0.0406 


0.9812 


4.1 


3.0 


5.7 


0.25 


0.3656 


0.8313 


2.32 


2 


03 


2 


52 


0.17 


0.0969 


0.9937 


3.25 


3.23 


3.2 


0.15 


0.3552 


0.9771 


2.39 


2 


38 


2 


24 


0.047 


0.1031 


0.9812 


3.09 


2.98 


3.3 


0.31 


0.3781 


0.9312 


2.34 


2 


31 


2 


24 


0.109 


0.1094 


0.9688 


2.98 


2.70 


3.44 


0.44 


0.4010 


0.8854 


2.31 


2 


23 


2 


26 


0.20 


0.1156 


0.9563 


2.91 


2.41 


3.62 


0.43 


0.4240 


0.8396 


2.30 


2 


15 


2 


31 


0.30 


0.1219 


0.9438 


2.86 


2.22 


3.70 


0.18 


0.4469 


0.7937 


2.27 


2 


06 


2 


36 


0.20 


0.1615 


0.9896 


2.80 


2.78 


2.7 


0.11 


0.4198 


0.9729 


2.35 


2 


34 


2 


19 


0.029 


0.1719 


0.9688 


2.70 


2.62 


2.77 


0.23 


0.4469 


0.9187 


2.30 


2 


28 


2 


17 


0.071 


0.1823 


0.9479 


2.63 


2.42 


2.87 


0.34 


0.4740 


0.8646 


2.28 


2 


23 


2 


18 


0.15 


0.1927 


0.9271 


2.58 


2.22 


3.00 


0.35 


0.5010 


0.8104 


2.27 


2 


17 


2 


20 


0.30 


0.2031 


0.9063 


2.54 


2.07 


3.07 


0.16 


0.5281 


0.7563 


2.25 


2 


11 


2 


26 


0.34 


0.2260 


0.9854 


2.59 


2.58 


2.48 


0.087 


0.4844 


0.9688 


2.32 


2 


32 


2 


17 


0.010 


0.2406 


0.9563 


2.52 


2.46 


2.50 


0.18 


0.5156 


0.9063 


2.28 


2 


27 


2 


14 


0.025 


0.2552 


0.9271 


2.46 


2.31 


2.58 


0.28 


0.5469 


0.8438 


2.27 


2 


24 


2 


13 


0.07 


0.2698 


0.8979 


2.43 


2.14 


2.67 


0.33 


0.5781 


0.7813 


2.25 


2 


21 


2 


14 


0.19 


0.2844 


0.8687 


2.41 


2.03 


2.74 


0.16 


0.6094 


0.7188 


2.24 


2 


16 


2 


18 


0.7 



7 



